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ABSTRACT 

We show that kicks generated by topological currents may be responsible for the large 
velocities seen in a number of pulsars. The majority of the kick builds up within the 
first second of the star's birth and generates a force about two orders of magnitude 
larger than a neutrino kick in the same temperature and magnetic field regime. Be- 
cause of the nature of the topological currents the star's cooling is not affected until 
it reaches 10 9 K; thereafter the current replaces neutrino emission as the dominant 
cooling process. A requirement for the kick to occur is a suitably thin crust on the 
star; this leads us to speculate that pulsars with large kicks are quark stars and those 
with small kicks are neutron stars. If true this would be an elegant way to distinguish 
quark stars from neutron stars. 
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1 INTRODUCTION 

The goal of this paper is to elaborate on a kick mecha- 



nism first discussed by Charbonneau & Zhitnitsky ( 2009 1 
that may explain pulsar velocities greater than 1000 km s -i . 
There have been a number of studies that have compiled and 
modelled the velocities of pulsars. Although they disagree 
on whether the distribution is indeed bimodal, they agree 
that a significant number of pulsars are travelling faster 
than can be attributed to neutrino kicks. The analysis of 
|Arzoumanian et al.| ( |2002| ) favours a bimodal velocity distri- 
bution with peaks at 90 kms -1 and 500 kms -1 with 15% 
of pulsars travelling at speeds greater than 1000 kms - . 
Alternatively |Hobbs et al.| ( |2005[ ) and |Faucher-Giguere fc| 
Kaspi ( 2006 1 both predict a single peaked distribution with 
an average velocity of ~ 400 kms -1 , but point out that 
the faster pulsars B2011+38 and B2224+64 have speeds of 
~ 1600 kms -1 . Large velocities are unambiguously con- 
firmed with the model independent mea surement of pul- 
sar B1508+55 moving at 1083±g° 3 kms -1 (IChatterjee et al 



2005) 



Currently no mechanism exists that can reliably kick 
the star hard enough to reach these velo cities. Asym met- 
ric explosions can only reach 200 kms -1 (Fryer 20041, and 
asymmetric neutrino emission is plagued by the problem 
that at temperatures high enough to produce the kick the 
neutrinos are trapped inside the star ( jSagert fc Schaffher-| 
Biclich 2008). Alterations of the neutrino model that take 



into account only a thin shell of neutrinos require large tem- 
peratures and huge surface magnetic fields. 
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The kick mechanism we will discuss relies on the exis- 
tence of topological vector currents of the form described by 



Charbonneau fc Zhitnitsky (20091, 
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where hr and are the chemical potential of the right 
and left-handed electrons, and $ is the magnetic flux. There 
are three requirements for topological vector currents to 
be present: an imbalance in left and right-handed particles 
A*l 7^ Hr, degenerate matter /i> T, and the presence of a 
background magnetic field B / 0. All of these are present 
in neutron and quark stars. The weak interaction, by which 
the star attains equilibrium, violates parity; particles cre- 
ated in this environment are primarily left-handed. The in- 
terior of the star is very dense, ~ 100 MeV, and cold, 
T ~ 1 MeV, such that the degeneracy condition /j, ^> T is 
met, and neutron stars are known to have huge surface mag- 
netic fields, B B ~ 10 12 G. Interested readers are encouraged 



to read Charbonneau & Zhitnitsky (2009) 



Kicks caused by the topological current are very similar 
to neutrino kicks. The helicity of the electrons and the sin- 
gle spin state in the lowest Landau level allow electrons to 
travel with a preferred direction through the neutron star. 
The difference between the two mechanisms is the linear 
momentum that the particle transfers to the star. Neutrinos 
have an energy close to the temperature of the star, while 
the electrons, which are degenerate, have an energy close to 
their Fermi momentum. We will discuss this critical differ- 
ence in Section [5] 

If the electrons carried by the current can transfer their 
momentum into space — either by being ejected or by radiat- 
ing photons — the current could push the star like a rocket. 
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In typical neutron stars this is unlikely because the crust 
(the region where fi ~ T) is thought to be about 1 km 
thick. Once it reaches this thick crust, the current will likely 
travel along it and be reabsorbed into the bulk of the star. 
But if the crust is very thin, the electrons may leave the 
system or emit photons that will carry their momentum to 
space. The crust for bare quark stars is thought to be about 
1000 fm ( Alcock et al.|1986} see also recent developments by 
Alford et al. (2006 1 ). With this in mind we conjecture that 
stars with very large kicks, v 3> 200 kms" 1 , are quark stars 
and that slow moving stars, v sC 200 kms -1 , are kicked by 
some other means, such as asymmetric explosions or neu- 
trino emission, and are typical neutron stars. Confirmation 
of this would provide an elegant way to discriminate between 
neutron stars and quark stars. 



2 SIMPLE MODEL OF QUARK STARS 

The kick will likely only occur in quark stars that are bare 
or have very thin crusts. We will consider only the existence 
of light quarks in the star, which attain equilibrium through 
the quark direct Urea processes, 



e + u (+ v e ) 
e~ + u (+ i> e ) 



d (+ v e ) 

S (+ V £ ) . 



(2) 
(3) 



where the neutrino terms only appear in the final state of 
the interaction. The direct Urea processes in quark matter 
do not have the same Fermi momentum restrictions as in 
neutron matter and are thus not suppressed. These processes 
are at equilibrium when, 



/J-u + Me = Pd and Pu + Pe — Ps 



(4) 



In this balance we ignore the neutrino chemical potential. In 
discussing the kick we will focus on temperatures where the 
star makes the transition from being opaque to neutrinos to 
being transparent. We will refer to this temperature range as 
the translucent regime. In the opaque regime the neutrino 
chemical potential is nonzero. As the star cools the mean 
free path of the neutrinos becomes large enough that all 
neutrinos created may freely leave the star. This is precisely 
the temperature at which the topological current begins to 
flow and creates a kick. In discussing these temperatures 
there is the concern that the neutrino chemical potential 
cannot be ignored. However, once in the transparent regime 
the neutrino chemical potential dissipates much faster than 
the time scale of the kick, thus we can safely ignore it. 
The quark matter must also be electrically neutral, 



2 2 1 

p/e = -n u - -rid - ^"s - n e = , (5) 

where p is the charge density and rii are the number densities 
of the quarks and electrons. We will consider the simplest 
model where the quark masses are set to zero, thus their 
Fermi momenta are equal to their Fermi energies. Unfortu- 
nately, this approximation requires the electron density to 
be zero to satisfy equation (jHJl, thus removing an essential 
element from the star: electrons. The strange quark mass is 
quite large and the number density should be described us- 
ing nonrelativistic Fermi statistics, which readmits elections 
into the mix to conserve charge. We will follow |Iwamoto| 
( 1982 1 by assuming that the quarks are all massless and an 



electron chemical potential is present to capture the con- 
tribution of the mass of the strange quark. The Fermi mo- 
mentum of the quark species can be written in terms of the 
baryon number density, denoted rib, as 
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MeV. 



(6) 



For typical densities in the core of the neutron star the 
Fermi momentum is k q ~ 400 MeV, which is greater than 
the strange quark mass. Neglecting the mass of the strange 
quark is a crude approximation but works well for this cal- 
culation. We can obtain the electron Fermi momentum us- 
ing the fraction of electrons to baryons, Y e = n e /nt, which 
yields 



1/3, 



73 



1/3 



MeV. 



Because fc e 3> m e , we find that k e ~ p e . The typical value 
for the electron fraction is Y e = 0.01. This value is not crit- 
ical in determining the strength of the kick because we will 
find that the dependance on k e cancels. However, as we will 
discuss in Section [5j the value of the Fermi momentum k e 
when compared to the core temperature of the star T cole is 
critical in understanding how topological currents generate 
large kicks. 



3 COOLING OF QUARK STARS 

In order to properly determine the size of a kick we need 
to understand how the core temperature of the quark star 
evolves with time. Unfortunately, kicks are likely to occur 
right after the birth of the star during the most poorly un- 
derstood stage of cooling. For insight into the cooling curve 
of a quark star we follow the work in two papers. The initial 



cooling of the star is described in Haensel et al. ( 1991 1, which 



focuses on neutrino diffusion through the star and thermal 
cooling. The star then cools until the neutrinos can escape 
the quark star and the cooling moves into a purely radiative 
regime as discussed in |Page fc Usov| ( |2002[ ). The part of the 
cooling curve between these two well defined mechanisms 
constitutes the translucent regime. 

The quark direct Urea processes, which are weak inter- 
actions, create both the neutrinos responsible for cooling the 
star and the parity that drives the current. The topological 
current only flows when the mean free path of the electron 
with respect to the weak interaction is larger than the neu- 
tron star. As they take part in the same processes, this is 
also when neutrinos start to escape the star. The mean free 
path of the weak interaction for an electron in quark mat- 



ter has been calculated in Charbonneau & Zhitnitsky ( 2009 \ 
and is 



1.2 x 10 1 



r c , 



10 9 K 



km . 



(8) 



Equation (JsJ tells us that for a star of radius R ~ 10 km the 
current starts to flow at T core ~ 7 x 10 10 K. At this point 
on the cooling curve the temperature changes on the scale 
of seconds. 

It is not clear how to model the short but critical cool- 
ing period that occurs in the translucent regime, as it is a 
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Figure 1. The dashed line indicates the translucent part of the 
cooling curve, modelled by exponential decay. The curve before 
the patch is take n from|Haensel et ah ( |1991[ ) and the curve after 
the patch is from [Page &: U sov (2002 . The black dot marks the 
start of the kick at t = 0. 



_L 



very complicated process, but we assume it must make a 
smooth transition from the primarily thermal cooling of the 



core found in Haensel et al. ( f 99f I to the neutrino cooling 
found in |Page fc Usov| ( |2002| ). Because we are interested 
in an estimate of the kick and not the exact details of the 
cooling we will interpolate between these two cooling curves 
by assuming the transition is smooth and exponential. We 
will use the cooling curve from Figure 2a of Haensel et al. 
( 1991 1 that models the cooling of the core of the star. We 



then piece it together using an exponential decay to cool the 
star to 10 10 ' 2 K and match the initialization of the cooling 



curve given in Figure 2 of Page & Usov ( 2002 1 . By requiring 
that the first derivatives are equal it is possible to find a 
unique exponential decay that patches these cooling curves 
together. Figure [I] plots this patched region that describes 
the translucent region of the cooling curve. 

It is possible that the cooling in the translucent regime 
is quite different from the exponential patch we propose. 
The severity of the exponential cooling provides us with a 
reasonable guess for the lower bound of the kick. Slower 
cooling in the translucent regime would result in larger kicks. 



4 ESTIMATING THE SIZE OF THE KICK 

The topological current runs along the magnetic field of the 
star and consists of electrons created in the core of the star 
travelling to the surface. Once at the surface these electrons 
transfer net momentum to the star. We assume that due to 
the thin crust of the quark star all of the electron's momen- 
tum is able to be transferred. The electrons may be physi- 
cally ejected from the star or the momentum may be trans- 
ferred by some radiative process such as bremsstrahlung. 



The size of the kick is directly related to the strength of 
the current, which is related to the creation rate of left- 
handed electrons in the quark star and the population of 
electrons in the lowest Landau level. The calculations re- 
garding the strength of the current in this paragraph are 



detailed in Charbonneau & Zhitnitsky (20091. We provide 



only a brief illustrative overview here. The direct Urea pro- 
cesses for quarks create electrons at a rate per unit volume 



w 

n 
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= 4.4 x 10 a 



10 9 K 



(9) 
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where a s — g 2 /4% « 0.4 and I ~ 20.6 is a dimensionless 
integral. The helicity of the electrons (A) = —0.84 and the 
fraction of electrons in the lowest Landau level are contained 
in the parameter, 
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The total number current for electrons reaching the surface 
of the star is given by, 
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(13) 



(14) 



where B c — 4.4 x 10 G is the critical magnetic field, T corc 
is the core temperature of the star, and no is nuclear density. 
The typical density for quark matter is ~ 10 no but could 
easily be higher. 

Because the current originates in the core of the star, 
we are interested in the strength of the internal magnetic 
field. The dipolar component o f the Vela pulsar at i ts surface 
is about B s ~ 3.5 x 10 12 G ( |Abney et aL]|l996[ ) but the 
field in the bulk of the star is likely much stronger based 



on virial theorem arguments in Lai & Shapiro (19911. The 



scalar virial theorem for the star states that 
2T + W + 311 + M = 0, 



(15) 



where T is the total rotational kinetic energy, LT is the in- 
ternal energy, W is the gravitational energy, and M is the 
magnetic energy. Because the kinetic energies are positive, 
the magnetic field energy can be at most equal to the gravi- 
tational energy, (47rB 3 /3)(S 2 nax /87r) ~ GM 2 /R, which yield 
possible core fields of _B max ~ 10 18 G. This is an extremely 
large field and is unlikely as it is a strict upper bound. 



Arguments by Abney et al. ( 1996 1 based on glitch timing 
of the Vela pulsar may put minimums on the core field 
depending on whether the protons are superconductors or 



not. As pointed out in Clem (19771, currents along type-II 



vortices are likely to destroy superconductivity. |CharborT| 
|neau fc Zhitnit sky (2007J also point out that the presence 
of topological currents could similarly destroy superconduc- 
tivity in neutron stars, thus the non-superfluid estimate for 
the internal field from |Abney et al.| |l996 ) is appropriate. 
In very dense matter they find minimum interior fields of 
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Figure 2. Time evolution of the kick for an internal magnetic 
field B = 10 B c . 



B n 



10 14 G. Based on this we choose a value of the core 



magnetic field to be Score = 10 B c 



The current, equation 1 13 1, gives a number rate of elec- 



trons reaching the surface of the star. To get the momentum 
transferred per second we simply multiply by the momen- 
tum each electron carries, k e , which is large due to the de- 
generacy and is given by equation Q. The acceleration for 
a IAMq star is then, 



^ e (^(-^eore)} 



1AM, 



© 



(16) 



We see that all the factors of fe e cancel because (J) contains 
a factor of k~ x . The kick is independent of the Fermi mo- 
mentum of the electrons, but as we will discuss later this 
does not mean the Fermi momentum is unimportant. 

As seen in Figure [2] the star quickly reaches a speed of 



1600 km s" 



(17) 



which is big enough to account for the large kicks seen in 
many pulsars. This result relies on the internal magnetic 
field of the star being two orders of magnitude larger than 
what we observe on the surface of most pulsars. We have 
also neglected the dissipation of electrons in the very thin 
crust of the quark star. As plotted, the entire kick seems to 
happen very quickly, but the current keeps running through- 
out the star's life. With a large internal magnetic field the 
mechanism can account for kicks seen in young pulsars such 
as Vela. But because the kick is constantly running, pulsars 
with smaller internal magnetic fields will eventually attain 
very large speeds very late in life. 



5 HOW ARE TOPOLOGICAL KICKS 

DIFFERENT THAN NEUTRINO KICKS? 

It is important to understand the differences between neu- 
trino kicks and topological kicks because superficially they 
seem very similar. In both cases a particle is created through 
the weak interaction with a favoured helicity (left-handedQ 
and the large magnetic field forces the electrons into Lan- 
dau levels. The lowest Landau level only accepts a spin 
down electron, whereas all the rest allow both spin up and 
spin down states. Combining the preferred helicity of the 
particles with the spin state asymmetry in the lowest Lan- 
dau level means that the neutrinos and electrons now travel 
with a preferred direction within the star. If the star is cool 
enough, these particles can escape the star without decay- 
ing. The neutrino does not participate in QED interactions 
and the electrons propagate due to the current, which is al- 
lowed because parity is a symmetry of QED. The particles 
that exit the star provide the linear momentum required to 
push the star. 

The electrons and neutrinos that contribute towards the 
kick are created at the same rate w, given by equation Q. 
We are only interested in interactions that create particles 
in the lowest Landau level so this rate is suppressed by a 
factor to account for the population of the lowest Landau 
level. The momentum transferred to provide the kick is sim- 
ply the momentum per particle that escapes the star. On 
the surface of the star, where the chemical potential is low, 
Boltzmann statistics are used to model the occupation of the 
lowest Landau level. This is what is seen in many neutrino 
kick models. This is because neutrinos require such hot tem- 
peratures to produce kicks that only a thin shell near the 
edge of the star is considered to contribute to the kick. In 
the core of the star, where the electrons are very degenerate, 
we must use Fermi statistics. This predicts an occupation of 



mi B 
~kfB~ c 



(18) 



in the lowest Landau level in the interior of the star. 

At core temperatures below 7 x 10 10 K the neutrinos 
can escape the star and the topological current is active. The 
flux of particles contributing to both the electron kick and 
the neutrino kick is about the same ~ nt, w. The difference 
between the two mechanisms comes from the momentum 
that the relevant particle carries. The neutrinos are created 
thermally and the typical momentum of a neutrino is equal 
to the temperature of the star T. The momentum of the 
electrons comes from the large chemical potential, /i e ~ 150 
MeV. The momentum transfer per unit time for neutrino^] 
is F v ~ Tn-^w and for electrons is F e ~ fif-n^w. When the 
kick starts the star has a temperature of only T ~ 1 MeV. 
The electron kick is stronger than the neutrino kick by a 
factor of 

h '"' (19) 



F v 



Me 



Initially, when the star is very hot, the electron kick is 

1 In a dense star the electrons are created with an average helicity 
(A) = —0.84, but for simplicity we will assume in this section that 
they are all left-handed. 

2 Notice energy per unit time has the T 6 dependance we expect 
for a calculation involving neutrino luminosity. 
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two orders of magnitude stronger than the neutrino kick. 
Furthermore, as the star cools the neutrino kick gets even 
weaker, while the electrons continue to have a momentum 
dictated by their chemical potential. This why electrons gen- 
erate larger kicks than neutrinos. 



6 HOW MUCH DOES THE CURRENT 

AFFECT THE COOLING OF THE STAR? 

At the beginning of the star's life the energy from the kick 
does not contribute to the cooling of the star, but later in life 
the current could become the dominant cooling mechanism. 
This may account for the colder than normal temperatures 
seen in the Vela pulsar. 

The reason is that only a small fraction of the electrons 
created in the neutron star actually escape due to the low 
occupation of the lowest Landau levels. As discussed earlier 
the estimate for the neutrino kick has a similar suppression 
factor, but all of the neutrinos actually escape the star, even 
those that do not contribute toward a neutrino kick. The 
electrons only propagate because the asymmetry in the low- 
est Landau level allows the helicity states to propagate out 
of the star. Those electrons that do not contribute toward 
the kick are trapped inside the star. The helicity states that 
reach the surface due to the current are the only electrons 
that contribute to the cooling of the star. The neutrinos cool 
the star with a luminosity L v ~ Tw where the electrons cool 
the star with an energy current (luminosity) of L e ~ k e rihW. 
The ratio of electron cooling to neutrino cooling is 



Page & Usov (20021. This patched cooling curve allowed us 



T 



At the star's birth the electrons cool the star at about 
1/100 the rate of neutrino cooling. As the star cools, even- 
tually L e /L v > 1 and the current becomes the dominant 
cooling mechanism of the star. This transition occurs at 
T t ~ 0.1 MeV ~ 10 9 K. The current is an additional cooling 
mechanism to consider in stars that have cooled below 10 9 
K. This may explain the unusually cold surface temperat ure 
of (7.85 ± 0.25) 10 5 K that is observed ( |Page et al.| 1996| ) in 
the Vela pulsar. 



7 SUMMARY 

We have found that the topological currents introduced in 



Charbonneau & Zhitnitsky (2009) can generate kicks greater 



than 1000 kms _T ? This may explain large pulsar kicks that 
are otherwise unattainable using traditional kick mecha- 
nisms. Topological kicks occur only if the electrons that 
make up the topological current can transfer their momen- 
tum through the star's crust. Because neutron stars have 
thick crusts, it is unlikely that they will be kicked due to this 
mechanism, but rather by mechanisms that generate weak 
kicks. In contrast, bare quark stars have very thin crusts 
that may allow the electrons to escape. Because of this dis- 
tinction we conjecture that many of the fastest pulsars are 
quark stars. 

In calculating the kick velocity we assumed that the 
cooling curve of the quarks star must make a smooth tran- 
sition from the primarily thermal cooling of the core found 



to integrate over time the momentum transferred from the 
topological current and obtain the strength of the kick. We 
have also demonstrated that the energy lost due to the cur- 
rent does not change the cooling of the star during the kick. 
But, after the star cools below 10 9 K the topological cur- 
rent becomes the dominant cooling mechanism rather than 
neutrino emission. 

The magnitude of the kick also depends on the strength 
of the magnetic field inside the star. The virial theorem puts 
a limit on the internal field of the star that is much higher 
than the surface field. To generate the kick we used a field 
strength of B COIC — 10 B c . This is a reasonable estimate for 
pulsars that have a surface magnetic field of B B ~ 10 12 G 
such as the Vela pulsar. 
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